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(<~^ \ It is generally assumed that the two dark components of the energy density of 

' the universe, a smooth component called dark energy and a fluid of nonrelativistic 

weakly interacting particles called dark matter, are independent of each other and 
Q i interact only through gravity. In this paper, we consider a class of models in which 

\ the dark matter and dark energy interact directly. The dark matter particle mass is 

■ proportional to the value of a scalar field, and the energy density of this scalar field 
comprises the dark energy. We study the phenomenology of these models and calculate 
the luminosity distance as a function of redshift and the CMB anisotropy spectrum for 

^ ' several cases. We find that the phenomenology of these models can differ significantly 

. from the standard case, and current observations can already rule out the simplest 

models. 
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1 Introduction 



Present observations of the universe strongly suggest that roughly ninety-six percent of the 
energy density of the universe is due to forms of matter and energy that are not described 
by the standard model of particle physics. These forms of matter and energy have little 
or no direct interaction with ordinary matter, and hence cannot be directly observed. The 
dark sector, called such because its energy and matter do not emit light, is typically divided 
into two categories, dark matter, which is clustered, and dark energy, which is smoothly 
distributed and presently causing the expansion of the universe to accelerate. It is also 
generally assumed that these two components of the dark sector are independent and do 
not interact directly. However, there are no experiments or observations that are known to 
explicitly preclude such an interaction. The goal of this work is to conduct a preliminary 
investigation into the constraints present observations place on possible interactions in the 
dark sector. 

Though there are presently no direct observations of the dark sector, there are neverthe- 
less many indirect measurements that give us clues about its nature. Measurements of the 
rotation curves of spiral galaxies, the temperature profiles of galaxy clusters, gravitational 
lensing of clusters, the large-scale motions of galaxies between clusters, and applying the 
virial theorem to clusters all require a mass for these objects much larger than that provided 
by the luminous matter |2l El IH El • It is widely believed that the extra mass is provided 
by a fluid of nonrelativistic, weakly interacting particles known as cold dark matter (CDM). 
Numerical studies of structure formation and comparison to statistical studies of galaxies 
and clusters support the CDM model on large scales [H]. 

Observations of the anisotropy spectrum of the cosmic microwave background radiation 
(CMB) jZlElE] and the magnitude-redshift relation for Type la supernovae (SNe la) [TUIITH 
IT^ IT^ indicate the presence of a second component of the dark sector, dark energy, which 
is smoothly distributed and presently causing the expansion of the universe to accelerate. 
The location of the first acoustic peak of the CMB anisotropy spectrum is predominantly 
governed by geometry of the universe, and hence the total energy density of the universe. The 
measured location of the first peak implies that the universe is very nearly flat, i.e. Euclidean. 
However, the measured clustered matter, including both ordinary matter and dark matter. 
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only comprises one-third of the energy density required for a flat universe implying the 
existence of a smoothly distributed dark energy. Measurements of the magnitude-redshift 
relation of SNe la indicate that the universe recently entered an era of accelerated expansion. 
The amount of dark energy required to cause the observed acceleration also makes up for 
the rest of the energy density needed to make the universe flat. 

The simplest explanation for the dark energy is a cosmological constant, A (for reviews 
see ^JEinni)- This hypothesis fits the data extremely well but faces significant theoretical 
questions. There is no known mechanism for producing a cosmological constant as small 
as is observed when compared with the Planck scale (A/Mj, ~ 10~^^°), and this constant 
must be very finely tuned to produce acceleration only very recently in the evolution of the 
universe. Another possible candidate is dynamical dark energy in which the dark energy is 
due to the potential energy of a scalar field, similar to the mechanism that is thought to drive 
inflation in the early universe [I71[Ii[inil2niEIlEai23Ell23l2Hll2Ill2Hl- A convenient 
parameterization of such models is the equation of state parameter, w = p/ p, the pressure 
divided by the energy density of the dark energy. For w ~ constant, the SNe la measurements 
and other observations restrict this parameter to the range —1.62 < w < —0.74 [21], which 
provides a stringent constraint on several interesting dynamical dark energy models. A 
cosmological constant would yield w = —1. If future observations indicate w < — 1, this 
would imply new physics in the dark sector, either a rapidly varying w [30J, a violation of 
a generally accepted condition on matter called the dominant energy condition [21], or a 
more complicated dark sector not yet explored. Moreover, though the CDM model with 
a cosmological constant, known as ACDM, impressively fits the current observations, new 
measurements in the near future will put the ACDM model through strict tests. It will be 
useful to have more general dark sector models around to both act as a foil to ACDM and 
as a possible replacement should it fail future tests. One possibility for a more complicated 
dark sector model is including an interaction between the dark matter and the dark energy. 

Casas, Garcia-Bellido, and Quiros considered models of scalar-tensor gravity in which 
the scalar coupled differently to different species of matter. Wetterich [^ pointed out that 
the scalar in such an interaction could constitute dark energy. Anderson and Carroll [S3] 
studied a cosmological model in which the dark matter particle mass was proportional to 
the expectation value of a scalar field. Similar models have also been studied by Bean [31] 
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and Farrar and Peebles jHSl- Amendola introduced a model with an exponential coupling 
between dark matter and a scalar field that acts as the dark energy [21]. This model has 
been extensively studied and shown to be consistent with present observations jSZl EHl IHH 
Hni E]- A particularly interesting property of this model is that it admits a late time 
attractor solution with Qdm/^de ~ 1/2, the ratio observed today. Other possible dark 
sector interactions have also been studied recently |121 UHl HH EEl EZl UHl HH] • Producing a 
consistent quantum theory of a dark matter - dark energy interaction may prove difficult j,45j, 
but in principle models could exist that can overcome these obstacles. Thus, understanding 
the phenomenology of such models continues to be of interest. 

In this paper, we consider the interaction studied by Anderson and Carroll [33 and 
recently by Farrar and Peebles The model discussed here has a direct coupling between 
the dark matter and dark energy. The dark matter particle mass is proportional to the value 
of a scalar field that acts as the dark energy. One should be concerned that by allowing 
the dark matter particle mass to change with time, the dark matter energy density will be 
too high or too low at earlier epochs. This could pose significant difficulties for structure 
formation and the CMB temperature anisotropy spectrum. The equations of motion for 
perturbations in these models differ, however from those in ACDM and other dynamical 
dark energy models. Thus, it is conceivable that there may be new effects in this model 
that compensate for the altered dark matter energy density. Careful calculations should be 
done to confirm or refute these concerns. The cosmological equations of motion governing 
this model for an arbitrary scalar field potential are derived in Sec. 2. In Sec. 3 we specify 
a potential and explore the phenomenology of the model, which can differ significantly from 
ACDM. We also calculate the luminosity distance relation for several models and compare 
the results to SNe la data. In Sec. 4 we consider perturbations in this model and calculate 
the CMB anisotropy spectrum for several cases. We conclude in Sec. 5 with a brief discussion 
of the results and future prospects for interacting dark sector models. 

2 Dark Sector Equations of Motion 

We begin this section with an intuitive motivation for the equations of motion governing a 
fiuid of dark matter particles with a rest mass that is proportional to the value of a scalar 
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field. A more detailed derivation from an action for two scalar fields follows. We consider 
here a model that is a flat, 3+1 dimensional Friedmann-Robertson- Walker (FRW) universe 
with metric 

^ _^^2 ^ a^Sijdx'dx^ . (1) 

The dark matter particles in this model will be similar to the particles of the CDM model: 
they will be collisionless and nonrelativisic. Hence, the pressure of this fluid vanishes and 
the energy density is given by the rest mass, m, multiplied by the number density, n, of the 
dark matter particles. We define 

m — A0, 

where A is a dimensionless constant and is a scalar field. The energy density and pressure 
associated with this fluid are thus 

Pdm — TTin — X4>n 
Pdm = 0. 

We will assume that this species of particle froze out in the early universe so that the 
comoving number density of dark matter particles is constant during the epochs of interest, 
i.e the particles are neither created nor destroyed. Thus, the number density is only a 
function of physical volume and n — noa~^ where no is the present number density of dark 
matter particles. 

One may be surprised that the dark matter is described by a pressureless fluid for which 
p 7^ po{a/ao)~^- The result p — po{a/ao)~^ for a pressureless fluid is due to the vanishing of 
the divergence of the stress-energy of a non-interacting pressureless fluid. The dark matter 
studied in this work interacts directly with a scalar field, and thus its stress-energy does not 
have vanishing divergence. We derive the evolution equation for the energy density of the 
interacting dark matter below. 

The scalar field that defines the mass of the dark matter particles has a standard kinetic 
term and a potential V{(f)). The energy density and pressure associated with the scalar field 
are 
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where an overdot represents a derivative with respect to time, d/dt. Since the energy den- 
sity of the dark matter particles depends on 0, the scalar field feels an additional effective 
potential when it is in a bath of dark matter particles. Taking this effect into account, the 
equation of motion is 

+ 3i/0 + ^ = O, (2) 

where H = a/a and 

V;ff(0) = V{(t>) + An0. 
Substituting this back into © yields 

+ 3i70 + ^ = -Anoa-3. (3) 
d(p 

The only difference between this equation and that for a noninteracting dynamical dark 
energy model is the term on the right hand side, which accounts for the interaction. 

In order to derive an evolution equation for the dark matter energy density, we first 
consider the divergence of the stress-energy tensor for each dark component. Since neither 
dark component interacts directly with any other species, the divergence of the sum of their 
stress-energy tensors must vanish. However, due to the interaction, the divergence of each 
stress-energy tensor is not necessarily zero. The derivative operator is linear, so 

which implies 

V^T(z)M)^ = -V^T(^)^. (4) 

The stress-energy tensor for the dark matter, T(^dm)'^u, is fairly simple. The only non- 
vanishing component is T(^dm)^o = —Rdm- For the scalar field, the stress-energy tensor 
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and its divergence is 



= - ( + 3if0 + ^ ) d,(j). (5) 
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Using the equation of motion for the scalar field Q, this expression simplifies to 

V^T(<^)^ = An9,0. (6) 

The evolution equation for the dark matter energy density is then calculated by combining 
dH) and ©: 

- ^fiTiDufo = PDM + 3-p = An0. (7) 

As a consistency check, we can use this equation and the definition of poM to determine the 
evolution of the number density as a function of the scale factor: 

Pdm = Ari0 + An0 + 3-An0 = An0, 

a 



which trivially reduces to 



n = —6—n 
a 



and is solved by n = ^0(0/00) ^5 as we had earlier assumed. 

The scale factor of the universe evolves according to the Friedmann equation 
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where p is the total energy density of the universe and G is Newton's constant. In the models 
considered in this paper, the universe contains the dark matter and dark energy discussed so 
far along with baryons, whose energy density scales as (a/ao)~^, and radiation, whose energy 
density scales as (a/ao)''^. Thus, the energy density as a function of the scale factor is: 

P = Pdm + P<p + PB,o(a/ao)~^ + PR,o{(^/ao)'^- 

We now continue with a more explicit derivation of these equations from an action for 
two interacting scalar fields, one of which will play the role of dark energy and the other 
dark matter. Consider a scalar dark matter particle, ip, and a homogeneous scalar field, 0, 
which serves as the source of the dark energy. The action for this model is 



S = d!^x 



(9) 



where R is the Ricci scalar associated with metric (0), is the determinant of the metric. 



and Mp = l/vovrG is the reduced Planck mass. Varying this action with respect to the 
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scalar fields yields 



+ 3if0 + ^ + A^c/)^^ = (10) 



- V^^ + 3ff7/^ + AVV = 0. (11) 

The spatial gradient term vanishes in the equation for because this field, which will make 
up the dark energy, is assumed to be spatially homogeneous. 

We now consider the equation of motion for the ip particle. The time scale associated 
with the dynamics of a particle of mass m is r ~ . To simplify the calculation, we will 
assume that ^ H ^ which is not a restrictive constraint during the epochs we will be 
considering. The mass of the ip field, = A0, evolves over short times as 

^ m^,* I 1 + ^^^^dt I 
- m,,(l + |-dt), 

where m^^^, is the ^/^ mass at some fixed time, 0=,, is the value of the field at the same fixed 
time, and dt is a time interval of order the time scale of the evolution of the -0 particle. We 
expect that the field, which is the dark energy in this model, evolves on a time scale that 
is comparable to the time scale associated with the evolution of the universe, i.e. 0/0 ~ if. 
Using this assumption and setting ~ r ~ yields 

m^^ IH . 

Since we have assumed that H/m^ <^ 1, the mass of the ip particle is a constant to a 
very good approximation. This reasoning also implies that the Hubble term in Eq. (fTTj) is 
negligible. This is reasonable because on time scales much shorter than that associated with 
the expansion of the universe, the expansion should be negligible. 

The equation of motion for the ip particle becomes after these approximations 

ip - V^%p + m^ip = 0, 

which not surprisingly is just the Minkowski space Klein-Gordon equation. The only differ- 
ence here is that over long times, i.e. on cosmological scales, changes. 



On timescales much smaller than the Hubble time, we can ignore the expansion of the 
universe and use the standard results from quantum field theory in Minkowski spacetime 
to study the properties of the particles associated with the ip field. We will work in the 
Schrodinger picture in which the field operators are time independent. We begin by decom- 
posing the ip field into its Fourier modes: 



where boldface indicates a 3-vector. The energy associated with this field is given by the 
associated (renormalized) Hamiltonian 



where Up = ylpP + and a|, and are the usual creation and annihilation operators. 
Since a^Op acting on a state yields the number of particles with momentum p, we can 
interpret (fT^ as 

E = J d^pNpUjp. 

Recalling that the dark matter particles are assumed to be nonrelativistic, i.e. the distribu- 
tion of ip particles is dominated by particles with |p| ^ m^, then it follows that the energy 
density in dark matter particles is given by 

Pi, = m^n^ = X4>n^. 

We now turn our attention to the equation of motion for the field, which we will 
treat classically. To do this, we need to interpret the \^(f)ip'^ term. Since we are treating the 
(f) field classically, we will consider the term: X^cf) {n\il)'^\n) where {n\il)'^\n) is the expectation 
value of in a state with n particles per unit volume. To evaluate the ip term, we begin by 
considering the ip'^ operator: 

i!p if!pL,i(p+p') >^ L (,. . J \ (,. ^ . J 



(27r)3 2ujp 



Since we are assuming that the ip particles are non-relativistic and it will greatly simplify 
the calculation, let us assume that all of the particles are in the p = state. Then we have 
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the (renormalized) result: 

{n\ip \n) = = = — 

Substituting this back into (fTUI) yields 

dV 

(P + 3H(P + — + Xn = 0, (13) 
a(p 

which is precisely the equation of motion motivated earlier. The rest of the classical discus- 
sion proceeds as already discussed in this section. Though we derived this equation for a 
scalar particle, it holds equally well for nonrelativistic fermions as shown in pB] . 

3 Cosmology with an inverse power law potential 

Exploring the phenomenology of the model described in the previous section requires a form 
for the potential of the scalar field, V{(f)). Since the mass of the dark matter particle depends 
linearly on 0, we choose a potential that blows up as approaches zero and thus prevents (p 
from becoming negative: 

V{(j)) = K(j)-''. (14) 

This form for the potential is also used in [32] and [33] . While this form of the potential is not 
necessary, the features of the cosmology resulting from this choice are significant, making it 
a good illustrative example. Substituting the potential (fTHl into the field equation Q yields 

+ 3i^0 - aK^-"-^ + Xnoa'^ = 0. (15) 

In the early universe, the effective potential looks like the upper left picture in Fig.^ The 
effective potential is very steep on both sides of the minimum, and the minimum is at a small 
value of 0. Due to the steepness of the effective potential and the effective friction due to the 
3H(j) term in ()15p. the field rapidly settles into the minimum of the effective potential for a 
wide range of initial conditions. As the universe expands, the number density of dark matter 
particles decreases, and the minimum of the effective potential moves to larger values of as 
seen in the upper right picture in Fig. ^ At the same time, the effective potential becomes 
shallower, and eventually the effective friction becomes important in the evolution of the 
scalar field. This slows the field so that it cannot keep up with the shifting minimum of the 
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effective potential. Eventually, the field is far from the minimum and effectively constant. 
This is shown in the lower right picture of Fig. Q The solid circle shows the actual field 
value, and the outlined circle shows the location of the minimum of the effective potential. 
When the field slows down sufficiently, the interacting dark matter particles begin to act just 
like ordinary dark matter particles with fixed mass and the field becomes like an ordinary 
non-interacting dynamical dark energy field. 

The evolution of the field when the field is in the minimum of the effective potential was 
first calculated in We reproduce and expand upon those results here. The value for the 
field at the minimum of the effective potential is obtained by solving for the field value that 
causes the first derivative of the effective potential to vanish, and the result is 

0mm = T — =00 — • (16) 

VAno/ \ao/ \ao/ 

In Fig. 121 we plot the value of the scalar field as a function of a/aQ for the case a = 1. In the 
early universe, the field takes the value 0min- At logiQ^a/aQ) ~ —0.6, the effective friction 
term begins to become important and slows the field down. This is seen in the plot as the 
actual value of the field falls below 0mm for logiQ^a/ao) > —0.6. 

While the field remains in the minimum of the effective potential, the mass of the dark 
matter particles is 

friDM = A -— — 



.XnoJ \aoJ 
and the energy density of the dark matter is 



PDM = Xno - — — . (17) 



We next consider the energy density associated with the scalar field. When the field is at 
the minimum of its effective potential, and hence following (|16p. its time derivative is 

= Y^H<p. 
1 + a 

The magnitude of the ratio of the kinetic to potential energy of the field is then 

V 4)-° ^ ' 
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Figure 1: The effective potential, Kff((/') changes as the universe expands due to the term 
that depends on the number density of dark matter particles. In all four pictures, the dotted 
line is the scalar field potential, the dashed line is the rest energy of the dark matter particles 
as a function of 0, and the solid line is the effective potential. In the early universe (top left), 
the number density of the dark matter particles is large causing the line associated with it to 
be steep. The field rapidly settles to the minimum of the effective potential. As the universe 
evolves (top right), this line becomes less steep, and the minimum moves to larger 0. The 
field (solid circle) follows the minimum until (lower left), the 3if0 term in the field equation 
(jS)) becomes comparable to the first derivative of the potential. After this time, the 3if0 
term acts as an effective friction and slows the field. Eventually the field is no longer in the 
minimum of the potential and is moving very slowly (lower right). The open circle in the 
lower right plot shows the minimum of the potential, and the solid circle shows the value of 
the field. At this point, the field is close to constant, and the interacting dark matter starts 
to behave as standard cold dark matter with a fixed mass. 
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Figure 2: Plot of log^Q^cp / Mp) versus logiQ^a/ao) for a = 1. The dashed hne shows the value 
of the field at the minimum of the effective potential, (pram/ Mp as a function of the scale 
factor. The solid line represents the actual evolution of the field. In the early universe, the 
field sits at the minimum of the effective potential, but eventually the effective friction due 
to the expansion of the universe causes the field to slow down and fall behind the shifting 
minimum of the effective potential. 
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Recall, however, that in order for the field to stay in the minimum of the effective potential, 
the first derivative of the effective potential must be much larger than the effective friction 
term in the equation of motion Q implying 

Comparing this to (fTH|) . we see that when the field is in the minimum of the effective potential, 
the kinetic energy is much smaller than the potential energy of the scalar field. Thus to a 
good approximation, the energy density in the scalar field is 

We see from fll7p and ((201) ^^at the dark matter energy density and the energy density of 
the scalar field have a constant ratio, 

Pdm 

= a, 

while the field is in the minimum of the effective potential. Substituting ()16|) into ()19|). we 
find that the approximation that the field is in the minimum of the effective potential breaks 
down when 

a* ( Ano\ 

Oo \aK ) 

An interesting parameter choice for this model is a = 1/2 and {Xn^/aK) ^ 1. In this 
case, there is an extended period of time during which pnul P4> = 1/2 as is observed. We 
plot the relative energy densities of radiation (Qr), baryons (VIb), dark matter (Qdm), and 
the scalar field (fi^) as a function of a/ao in the bottom panel of Fig. El In the top panel, 
we plot the relative energy densities for a ACDM model for comparison. 

For any other value of a, we require a*/ao < 1 in order to have Pdm/ P<t> ~ 2 today. In 
Fig. m we plot the relative energy densities in a number of models with different values of a 
and fixed Xtiq. Changing Xtiq does not significantly alter these plots, since we must adjust 
the value of K in order to get the right energy density in baryons and radiation today. 

It is typical in models with small values of a that the universe is baryon dominated 
during structure formation, which would likely produce structure that is significantly different 
from what we observe. Calculating the formation of structure in these interacting dark 
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Figure 3: The bottom panel is a plot of the relative energy density in radiation {Qr, long- 
dashed line), baryons {fls, short-dashed line), dark matter (fi^jM, dotted line), and the scalar 
field acting as the dark energy {Q^, solid line) for a = 1/2 and (Xno/aK) ^ 1. The top 
panel is the same plot for ACDM with Q^, the relative energy density in the cosmo logical 
constant, replacing and is included for comparison. Eventually, the dark energy will 
dominate in the interacting model, but this can be put off as long as desired by increasing 
(Xno/aK). Note that in this case, Qe,m/^<i> = 1/2 is sustained for many e-foldings of the 
scale factor. Another feature of this model is that the onset of dark energy density takes 
place over a much longer time than in the ACDM scenario. Finally, the universe is never 
dominated by dark matter, but instead goes through an epoch of baryon domination. 
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a=l/3 a=Z 




-4 -2 2 -4 -2 2 

logio(a/a„) log((,(a/a„) 



Figure 4: Plots of the relative energy density in radiation long-dashed line), baryons 
{VLb, short-dashed line), dark matter {VLdm^ dotted line), and the scalar field acting as the 
dark energy {VL^, solid line) vs. logiQ^a/ao) for various values of a. 
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sector models is beyond the scope of this paper. Instead, we will compare these models to 
observations in two other areas, the magnitude-redshift relation of SNe la in this section and 
the CMB anisotropy spectrum in the next. 

Often in comparing dynamical dark energy models to SNe la data, one calculates the dark 
energy equation of state parameter, and compares this theoretical value to the values for 
w = constant dynamical dark energy models allowed by the data. Such a comparison will not 
work in this case, however, because the w = constant dynamical dark energy models assume 
that the dark matter energy density scales as (a/ao)~^, which is not the case in interacting 
dark sector models. Thus, we must calculate the predicted magnitude-redshift relation for 
SNe la for each model and compare the predicted relation to the measured one. 

The magnitude-redshift relation of SNe la measures the luminosity distance, di, as a 
function of redshift, z, via the relation 



where m is the relative magnitude of the supernova and M is its absolute magnitude. Typ- 
ically, the results of these measurements are plotted as the difference between the observed 
(m — M) and that expected from a standard model called the "empty" universe: 



The luminosity distance is computed from the evolution of the scale factor through the 
relation 



We have calculated the theoretical A(m — M) for a number of models, and the results are 
plotted in Fig. El along with binned data. For comparison, we have also plotted the ACDM 
result. Note that large values of a do not fit the data well because they do not have enough 
acceleration in the expansion of the universe to match the data. Small values of a may 
also be problematic because acceleration begins too early to match the result of SN1997ff 



m- M = logio dh^z) 



A(m - M) = logio ddz) - \og,, dr^''{z). 




and the result for the "empty" universe is 




z{l + z/2) 
Ho 
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at z ~ 1.7. However, further high-redshift data will be required to strictly rule out such 
models. Similar results have also been obtained for the models considered by Amendola and 
his group ISO]. 



4 Perturbations and the CMB anisotropy spectrum 

In this section we derive the equations of motion for perturbations of the scalar field and the 
dark matter particles, which will be used in calculating the CMB anisotropy spectrum. We 
will work in conformal time [adr = dt) and in the synchronous gauge where the metric is 



The scalar field and the dark matter energy density can be written as the sum of a back- 
ground, average value and a perturbation that is assumed to be small when compared with 
the background value: 



PDM — PDM + SpDM- 

The dark matter perturbation has two components, the perturbation in the number density 
of the dark matter particles, Sn, and the perturbation in the scalar field, 6(f): 

8pDM = A0(5n + Xn6(f), 

where n is the average number density of dark matter particles. It will be useful to also 
define the relative perturbation in the number density of the dark matter particles: 

5n 

On = —■ 

n 

The equations of motion for the background are d?)), and (jH}, which we reproduce 
here in terms of the conformal time, r: 



ds'^ = a^(r) -dr'^ + {5ij + hij) dx'dx^ 



= + (5</) 



<p + 2—(p + a —J = —Xrio/a 




(21) 



p'dm + ^—PDM = An0' 



(22) 



a 
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Figure 5: Comparison of several models with SNe la data. The solid lines lines represent 
models with (from top to bottom) a = 1/2,1,2,5,8. The theoretical curve for ACDM is 
shown as a dotted line for reference. The solid circles are binned SNe la data found in |12j, 
and the point at z = 1.7 is the upper limit for SN1997ff from The models with a > 2 
are disfavored by moderate redshift observations. 
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- 

where a prime indicates a derivative with respect to the conformal time, d/dr. The total 
scalar field satisfies the equation 

dV 

+ — + An = 0. 

d(p 

A Fourier mode of the scalar field perturbation 



/2n 

satisfies the equation of motion 

6(j)'l + 2-50', + (e + a^^] + Xna^dn = (24) 
a \ dcp'^ J 1 

where we have used the background equation of motion (PT|) and h is the trace of the metric 
perturbation hiy From here on we will suppress the subscript k on the perturbed variables. 

In order to get the equation of motion for the dark matter perturbation, we will need to 
consider the divergence of the perturbed stress-energy tensor. Let us begin by considering 
the perturbed stress-energy tensor for the scalar field, 

5r(^)^ = {dn<\>) (9,0) + (9^0) {dM) - 5\ Q {dJ<P) + ^ (dj) (9"50) + ^50^ , 

and its divergence, 

W^5T^^f, = 1^50" + 2^50' + kHcj) + + \h'^'^ 9,0 

\ a d(p J 

= An5„ (9,0) + An (9,50) , (25) 

where in the last line we used (PT|) and (j21I). 

The perturbations to the dark matter stress-energy tensor are defined to be (following 

EH) 

STdm^O = —SpDM 
STdm^i = —STdm^O = pDMVi 
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The other components of the perturbed stress-energy tensor vanish. The divergence of the 
stress-energy of the combined dark sector must vanish to all orders. Given (j2Hl), "we demand 



The resulting equations of motion for the dark matter perturbations are 

K = -<>-\h' (26) 

, . .hl,i\e,%'^, (27) 

\a (p J a"^ (p 

where 6 = ik^vj. 

Equations ((211), and (j77|l were integrated along with the standard equations for 
baryons and photons using a modified version of CMBFAST ^2|- The results are plotted in 
Fig. El We were unable to fit any of the models to the WMAP CMB anisotropy spectrum [H31 
with values of the cosmological parameters that are consistent with other experiments. In the 
first five panels of Fig. El we show for illustrative purposes a theoretical CMB temperature 
anisotropy spectrum with the cosmological parameters set to the best fit values in jHlj except 
for allowing a variation in the scalar spectral index, Ug, a parameter that describes the initial 
density power spectrum. Cosmic variance is shown as dotted lines, and the error bars on the 
binned points include only experimental error. Effective values were calculated for these 
models following [53] , and were all many orders of magnitude larger than the best fit ACDM 
model. In the lower right panel, we show all five spectra evaluated with = 1 in order to 
illustrate the dependence of the spectrum on a. 

Much of the error in these plots is due to discrepancies between the theory and the data 
around the first peak in the spectrum. One can imagine, though, that a better fit around the 
first peak could be obtained by performing a fit over a larger number of parameters. Detailed 
parameter fitting, though, is beyond the scope of this paper. Even with a better fit around 
the first peak in the spectrum, these models face two qualitative problems in predicting the 
CMB temperature anisotropy spectrum, one at small scales (large /) due to a low value for 
the dark matter energy density at last scattering and the other at large scales (small /) due 
to an enhanced ISW effect (for a review of CMB physics, see |[55|). 

It is not surprising that these spectra look like those for a low dark matter density 
universe. In the models considered here, the dark matter particle mass is increasing with 
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Figure 6: Plots of the CMB temperature anisotropy spectrum for several models. In the 
first five panels, the theoretical spectrum is compared to the WMAP results [HS]- For each 
model, the best value of the scalar spectral index, n^, is shown. The dotted lines represent 
cosmic variance. The other cosmological parameters, are set to the best fit values in for 
illustrative purposes. In the lower right panel, we plot all five spectra with the same Us = 1, 
though for a = 1/2, 1, 2, the differences between the spectra are less than the line thickness 
in this plot. The solid line represents these values of a, the long dashed line is a = 5, and 
the dot-dashed line is a = 8. 
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time. Holding the present dark matter energy density constant, a larger change in the dark 
matter particle mass since last scattering implies a lower value for the dark matter energy 
density at last scattering. At small scales, the spectra for these interacting models look very 
much like the spectra for ACDM models with a low value for Qdm, models which also have 
a low dark matter energy density at last scattering. 

The ISW effect, which affects the power in the anisotropy spectrum at large scales, is 
due to the integrated evolution of the gravitational potential along the path of the photons 
from the last scattering surface. During a matter dominated phase, large scale density 
perturbations evolve in a manner that keeps the gravitational potential constant. Thus, in 
a universe that is matter dominated from the time of last scattering to the present, no ISW 
effect is observed. During periods when the dominant equation of state parameter, the ratio 
of the pressure to the energy density of the dominant energy component of the universe Wd, 
is changing, the large scale gravitational potential also changes 'STj. Thus, the ISW effect 
is observed in the ACDM model and in dynamical dark energy models due to the recent 
onset of vacuum domination and the resulting change in Wd from zero to —1 58 . In the 
interacting dark sector models discussed here, the onset of vacuum domination takes place 
over a long time, and hence Wd is also changing over an extended time producing a large 
ISW effect. Moreover, the mass of the dark matter particles is increasing with time, which 
also causes the gravitational potential to change, and hence may cause an increase in the 
observed ISW effect. 

The dependence of the anisotropy spectrum on a shown in the lower right panel of Fig. IHl 
is predominantly due to the dependence of the ISW effect on a. When the scalar field is 
in the minimum of its effective potential, the energy density in both dark matter and dark 
energy scale as (a/ao)"^"/''^^"^ This is equivalent to a model with a single fluid with an 
energy density evolving as (a / aQ)~^^^~^'^'> with w = —1/(1 + a). For larger values of a, w is 
closer to zero, the value of w^s during baryon domination. Hence, larger values of a imply a 
smaller change in Weff as the universe transitions from a baryon dominated epoch to a dark 
matter/dark energy dominated epoch, which in turn implies a smaller ISW effect. In the 
spectra shown in Fig. El the ISW effect is larger for the smaller values of a. 
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5 Conclusions 



The main result of this work is a demonstration of two challenges facing cosmological mod- 
els with an interaction between the dark matter and dark energy. At the outset we were 
concerned that an interaction that causes the mass of the dark matter particles to change 
in time would lead to difficulties in reproducing the measured temperature anisotropy spec- 
trum in the CMB. In the models considered here, the mass is a monotonically increasing 
function of time. As a result, the energy density of the dark matter diminishes more slowly 
in this model than in the ACDM model, which in turn implies a lower ratio of dark matter 
energy density to radiation energy density at the time of last scattering than in ACDM. 
It was conceivable, though, that since the equations of motion for the perturbations in an 
interacting model differ from those in ACDM, some new effect could compensate for the low 
dark matter to radiation ratio. Having performed the calculation, we now know that this is 
not the case, and the initial assumption that the low dark matter to radiation ratio would 
be seen in the CMB temperature anisotropy spectrum was correct. 

It is not difficult, however, to imagine models that can overcome this challenge. One 
possibility is a model in which the dark matter particle mass at the time of last scattering is 
comparable to the dark matter particle mass today, but varied during the intervening time. 
For example the dark matter particle mass could have been decreasing from a large value in 
the early universe, reached a minimum at some time between last scattering and today, and 
be increasing today. Such a model could not be ruled out by the considerations in this work, 
and detailed study of structure formation in such a model would be necessary to test it. 

Another possibility is to finely tune the initial conditions so that the scalar field does not 
evolve enough from the time of last scattering to today to be ruled out by the CMB data. 
In the models considered here, generic initial conditions lead to the scalar field sitting in the 
minimum of its effective potential for a significant time. One could finely tune the initial 
conditions so that the scalar field does not reach its minimum, and the effective friction 
due to the expansion of the universe prevents the field from evolving too much so that it 
can satisfy the CMB constraints. Such a situation, though aesthetically displeasing, still 
remains a possibility. A similar effect can also be obtained by either introducing a new dark 
sector field or a different potential that holds the field near a fixed value. Ruling out these 
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situations would again require a more detailed study of structure formation, and this is the 
approach of |35|. 

The second challenge to interacting models is avoiding an enhanced ISW effect at small / 
in the CMB temperature anisotropy spectrum. One of the potential aesthetic benefits of an 
interacting dark sector model is that the transition from a dark matter dominated epoch to 
a vacuum dominated epoch takes place over an extended period of time when compared to 
ACDM, thus softening the coincidence problem. Unfortunately, this work has demonstrated 
that, at least in the specific model considered here, allowing the dominant effective equation 
of state parameter to vary from zero during matter domination to nearly —1 during vacuum 
domination over too long a time leads to an enhanced ISW effect that is in conflict with 
observation. The methods discussed above to overcome the challenge of a time varying dark 
matter particle mass would likely also overcome this challenge as well. Such models allow for 
only small deviations from ACDM and more standard dynamical dark energy models, and 
hence the evolution of the effective equation of state parameter would be similar to ACDM 
and not differ from what is observed. 

In this work we have used the CMB temperature anisotropy spectrum to rule out a class 
of cosmological models with an interaction between the dark matter and dark energy. These 
models, however, are only among the simplest possibilities for such an interaction. More 
complicated models may well overcome the challenges pointed out here, and with strict new 
observational tests of the ACDM model on the way, it will be useful to continue to study 
interacting dark sector models. 
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